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We provide a systematic treatment of Ai-photon transport in a waveguide coupled to a local 
system, using the input-output formalism. The main result of the paper is a general connection 
between the Wphoton S matrix and the Green functions of the local system. We also show that the 
computation can be significantly simplified, by exploiting the connectedness structure of both the S 
matrix and the Green function, and by computing the Green function using an effective Hamiltonian 
that involves only the degrees of freedom of the local system. We illustrate our formalism by 
computing A^-photon transport through a cavity containing a medium with Kerr nonlinearity, with 
N up to 3. 

PACS numbers: 


I. INTRODUCTION 

The capability to create strong photon-photon interaction at a few-photon level in integrated photonic systems is 
of central importance for quantum information processing. To achieve such a capability, an important approach is to 
use the so-called waveguide quantum electrodynamics (QED) system, where one confines the photons to a waveguide 
that is strongly coupled to a local quantum system. Experimentally, the waveguides that have been used for this 
purpose include optical fibers [1], metallic plasmonic nanowires [2], photonic crystal waveguides [3], and microwave 
transmission line [3] . The local quantum system typically incorporates a variety of quantum multi-level systems such 
as actual atoms quantum dots [2 E], or microwave qubits @], where the strong nonlinearity of these multi-level 
systems forms the basis for strong photon-photon interactions. These multi-level systems moreover can be embedded 
in cavity structures to further control their nonlinear properties [SHin]- 

The rapid experimental developments, in turn, have motivated significant theoretical efforts. From a fundamental 
physics perspective, the photon-photon interaction is characterized by the multi-photon scattering matrix (S matrix). 
Therefore, a natural objective for theoretical works is to compute such multi-photon S matrix. Moreover, from 
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an engineering perspective, the systems considered here are envisioned as devices that process quantum states. To 
describe these systems as a device one naturally have to specify its input-output relation. The S matrix, which relates 
the input and output states, therefore provides a natural basis for device engineering as well. 

Motivated by both the physics and engineering considerations as discussed above, a large body of theoretical works 
have been therefore devoted to the computation the S matrix of various waveguide QED systems [TTHSl]. These 
computations, however, are limited in two important aspects: 

(1) All of these computations are carried out for a specific local quantum system. In most of these cases, the 
methods that were used were tailored to the property of the specific system. In the wavefunction approach for S 
matrix calculation [iiiiiiiiiiisiiiHiiiiiiisiiisiiiHiisniEa, for example, the ansatz for the wavefunction used is 
specifically related to the local quantum system. As a result, it has been difficult, from these calculations, to identify 
the general features of S matrices for waveguide QED systems. 

(2) With a few exceptions [TH ITH] . almost all previous calculations have been carried out for either single or two- 
photon S matrix. On the other hand, in quantum information processing, there is a strong effort to create and 
understand highly entangled states with more than two photons |33H35j . It is important to understand whether 
waveguide QED system can be used for such a purpose. Thus, computation of Wphoton S matrix with iV > 2 is 
essential. 

In this paper, we extend the input-output formalism [m HOI im 1311 [35] to provide a systematic computation of 
Wphoton S matrix for waveguide QED system. The main result of the paper is the relation between the 7V-photon 
scattering matrix, and the Green function of the local system. We prove this result using only a quantum causality 
condition, without the need of knowing the specific details of the Hamiltonian of the local system. The main result 
is therefore generally applicable for a large number of waveguide QED systems with different local quantum systems. 
We also discuss the general connectedness structure of both the S matrix and the Green function, which arise from 
the local nature of the interaction, and show that such connectedness structure can be used to significantly simplify 
the computation. 

Our work represents a significant step forward in the understanding of waveguide QED system. Our results here 
highlight some of the universal nature of the properties of these strongly correlated systems that has not been 
emphasized before. As a computational method, our work here leads to an approach for systematic computation for 
Wphoton scattering matrix that is directly applicable to a large number of different systems. 
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Our work utilizes the input-output formalism developed in standard quantum optics literature. However, our focus 
here is different. Whereas much of the standard quantum optics literature have focused on computing properties 
related to an input state that is a coherent state, a thermal state, or a squeezed state, here we focus exclusively on 
computations for Fock state input. In general, the transport property of Fock states is qualitatively different from that 
of the coherent state. As a prominent recent example as developed in the context of Boson sampling problem [37113H], 
it has now been recognized that the A^-photon Fock state transport in a linear waveguide network is computationally 
hard |39j . even though the transport properties of the same network for coherent states are well known. Similarly, 
in the system that we are discussing here, while many properties of the system in the presence of a coherent state 
input can be and have been computed with standard quantum optics tools, much less is known about how to compute 
transport properties of the same system with a Wphoton Fock state input. 

The paper is organized as follows. In Section II we briefly review the input-output formalism and derive the 
quantum causality condition. In Section III we prove certain time-ordering relations, which are the key to compute 
the S matrix of waveguide photons. In Section IV we prove the connection between the V-photon S matrix and 
the time-ordered local system’s Green function. This derivation represents the main result of the paper. To further 
simplify the calculation, we study the connectedness structure of S matrix in Section V. We also show in Section VI 
that the system’s Green function can be computed with an effective Hamiltonian approach. Finally, in Section VH 
as an example of the application of this formalism, we calculate the exact V-photon S matrix up to V = 3, when the 
local system is a cavity containing a medium with Kerr nonlinearity. 


II. A BRIEF REVIEW OF THE INPUT-OUTPUT FORMALISM 


We start with a brief review of the input-output formalism, highlighting only those aspects that will be required for 
the paper here. More details can be found in [171136]. Following [17[[36], we consider the Hamiltonian of a one-mode 
waveguide coupled to a local system with finite degrees of freedom {h= 1): 


H = 



( 1 ) 


where ^ is the coupling constant between the waveguide and the system, and is assumed to be frequency independent. 
Cfc (c^) is the annihilation (creation) operator of the photon state in the waveguide satisfying the standard commutation 
relation [ck, c^,] = 5{k—k'). We consider only a narrow range of frequencies, in which the waveguide dispersion relation 
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can be linearized, and the group velocity of the waveguide is taken to be 1 . a(a^) is one of several possible system 
operators that are assumed to commute with Ck, c|,. In this section we assume a to be arbitrary. In practice a can 
either be a bosonic operator describing a cavity mode [n m [551 nHi ng, or a spin operator for atom-waveguide 
interaction [HdiiiiniiMiEni- 

The aim of the paper is to develop a systematic approach to compute the IV-photon scattering matrix. In general, 
the A^-photon S matrix is related to the input and output operators by m 


^ ^ f dt 


- \J[ J ^ 

\ 2=1 '' J = 1 


3 g 


N N 

i=i j=i 


( 2 ) 


in which we define the input and output operators as: 


Cin(i) = J dkck{to)i 

Cout(^) — dk Ck(tl^ ^ 




ifc(t-ti) 


(3) 


with to “t —oo ,ti —>■ - 1 - 00 . We note that Cin(t) and Cout(i) consist of Heisenberg operators of waveguide photons at 
time —oo and -l-oo, respectively. They also satisfy the commutation relations 


[cin(t),Ci„(t')] = [Cout(i))Cout(i')] = 0: Ci„(t),c|'„(t') = Cout (t), (^') = d{t - t'). 


4 - 


(4) 


Following the standard procedure [T71|35], one can develop the standard input-output formalism that relates Cin, 
Cout and a, as: 


Cout(0 = Cin(t) - *V7a(t), 


(5) 


^ V7 [a. Cin - i 


u •T' t ' 

a, Hsys - i-a'a 


( 6 ) 


or 


da 


— = -i [a, a'l'] Cout - i 


a, Hsys + i-a'a 


where 7 = 27r^^. Integrating § and Q from t = —00 and t = 00 , respectively, result in: 


(7) 


i{t) = a(— 00 ) — i dr 


i{t) = a(-l-oo) — i dr 

J + 00 


u ■'>' t ■ 
a, Hsys - i-a'a 


u , ■7' t ' 
a, Hsys + I—a' a 


-iy/j dr [ 0 , 0 '^'] Cin 


- * Vt / dr [a, a'l’] Cout ■ 
J + 00 


( 8 ) 

( 9 ) 
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where the integrands in ([^ and are operators at time r. 

(|^ and § can be used to prove a quantum causality relation. When using ^ to evaluate a{t) or a^{t), the integral 
in should result in an expression that involves only Cin('r) and cI^(t) with t < t. Therefore, by the commutation 
relation Q above, one concludes from Q that for t 

[a(t), m] = [a(-oo), , [a\t), I{t')] = [at(-oo), , (10) 

where I{t') is a shorthand notation for the input operators that represent either Cin(t') or On the other hand, 

the operator I is really a Heisenberg operator at time —oo as can be seen in ([^ above, and hence commute with 
a(—oo) and a^(—oo). Therefore, we have 

[a(t),/(t')] = = 0 : for t ^ t', (11) 


Similarly, one can prove 


[a(t), 0{t')\ = [a^(t), 0{t')\ = 0, for t ^ t', 


( 12 ) 


where 0{t') is a shorthand notation for the output operators that represent either Cout(t0 ^^^^(t'), by utilizing ([9 
and the fact that the output operators are really Heisenberg operators at time +oo. Following [36j, we refer to 0 


and (12) as the quantum causality condition. The operator a(t), which characterizes the physical field in the local 


system, depends only on the input field Cin('r) with r < t, and generate only output field Cout(T) with t > t. 


Ill RELATION INVOLVED TIME-ORDERED PRODUCT 


Having reviewed some of the basic aspects of the input-output formalism, we now proceed to compute the N- 
photon S matrix as defined in ([^. For this purpose, we first consider some of the properties of a time-ordered product 
involving a and the input or output operators. We note that: 


Ta{t)I{t') 

= a(t)/(t'), 

(13) 

Ta{t)0{t') 

= 0{t')a{t). 

(14) 


Take (13) as an example, by definition, Ta(t)I(t') 


a{t)I{t') for t ^ t'; When t < t', by the quantum causality 


condition of (11), we have Ta{t)I{t') = I{t')a{t) = a{t)I{t'), completing the proof. (14) can be proved similarly. 
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More generally, we have the following relation regarding the time-ordered product: 


'iJ 


TY[a{U) 

i 

rno(i') 


rllaf' 

3 

T\{a{U) 


(15) 

(16) 


where bracket is used to indicate the range over which the time-ordering is being applied. 


(15) and (16) can be proved in a similar way. Here we show only the proof of (15). The proof of (15) can be 


constructed from induction with respect to the number of operators. The base case is already proved in (13). Now 


suppose (15) holds for all cases involving a total number of N operators of a and I. Consider a time-ordered product 


involving iV -|- 1 operators, if the operator with the largest time label is a{t„ 


T ]^a(f/)/(<') = a{tn 

1,3 






— Q-(tmax) 


TX{a{U) 




r J|a(t/) 




(17) 


where the definition of time-ordered product is used in the first and last steps and the induction hypothesis is used 
in the second step. On the other hand, if the operator with the largest time label is /(tmax); we have 


i,J 




rY[a{t 


lit max) 


— / (^ma: 


T a(t 

i 

T\{a{U) 




r\{iit'j) 

j 


(18) 


where we use the induction hypothesis in the second step and the commutation relation 0 in the third step. 


Therefore, (15) holds for iV -|- 1 operators, completing the proof. 


IV JV-PHOTON S MATRIX 


Using (15) and (16), we now compute the 7V-photon S matrix defined in (2). We first evaluate its Fourier transfor¬ 
mation in the time domain 

N N 

= (0| Cout(ii) c|,,(tj)|0) . 
i=l j=l 


( 19 ) 
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From now on, we refer both ([^ and (19) as S matrix. The computation is as follows: 


5, 




N 


N 


J i=i 


Cout(ti) 

. 2^1 
N 

TY[{Cinit’i) - iy/j 


2=1 


N 




i=i 


N 


= 

Ivl—Q 
N 


M 


'rn«( 4 „(«)) 


S = 1 

M 


N-M 


N 




M 


i-1 

N-M 


T Cin(ts^(„)) 

m—1 
t 

M—0 Bm,Dm s —1 2"=1 P m—1 

N MM 

M—0 Bm,Dm s —1 1 B m—1 

N MM N-M 

= E<-^)" E (oi^n“(*'i.«(.))n“'fe«w)io)i: n (20) 




N-M 


M=0 


Ba4 ,Da 


s=l 


r—1 


P m—1 


In the first step above, we introduce the time-ordering operation for the product of Cout operators since these 
operators commute with one another. 

In the second step, we use the input-output relation ©■ 

In the third step, we expand the product in the bracket. In this derivation, for a given subset B of {l...A'^}, we use 
B‘^ to represent its corresponding complementary subset, and B{s) to represent its s-th element. The summation is 


over all subsets of {1...7V}. Bm is a subset with M elements. We use the time-ordering relation (15) to place all Cm 
operators to the right of the a operators. 

In the fourth step, we first remove the time-ordering operation for the product of Cjn and then contract the Cjn 
operators with the cl^ operators. Since there are N — M Cin and N cl^ operators, in each term in the resulting 
summation we select M operators as indexed by a M-element subset Dm, and perform a full contraction using 
the remaining N — M operators. Such contraction results in a summation over all possible permutations P of D'^. 


Finally, we restore the time-ordering operation on the product of M operators. The use of (15) then results in a 
time-ordered product of all the a and c-^ operators. 

In the fifth step, we again use the input-output relation ([^. 


In the last step, we expand the product involving a^, and then apply (16) to each term in the product expansion. 


For every term that contains at least one operators, the use of (|16[) resulting in such output operators being 


placed on the left-most positions of the operator product, and hence such a term vanishes. 















Therefore, we obtain the first main result of the paper: 


N 


M 


M 


N-M 


= E (oi'rn«(4„(.))n“'W))io)E n 




( 21 ) 


M—0 Bm,Dm s —1 

We define the time-ordered 2M-point Green function 


P m—1 


Git'i ■ = {0\Ta(t[) ■ ■ ■ a(t'J^^)a\tl) ■ ■■a'^itM)\0) 


( 22 ) 


and its Fourier transformation 


M 


G(pi • • -pM', fcl ■ • • kM) = P[ 


dt'i 


M 

P[ 


dti 


1 e I G{t'i ■ ■ ■ ti ■ ■ ■ Im) ■ 




(23) 


From (211, the A^-photon S matrix ([^ in the frequency domain is 

N 

E E ^ iPBMTkDM)'^2i IT ^ (pBuim) — kpDl^{Tn)j j 


N-M 


S; 


'pi---pN',ki---kN — 


(24) 


M—0 BmiDm 


p m—1 


where we use the shorthand notation pbm = {Pi|* G Bm} and = {ki\i G Dm}- In (24), represents the 


frequencies of the incoming photons that bypass the local system. These photons do not change their frequencies, as 
signified by the (5 functions in ( |24[ ). Whereas kp^ represents the frequencies of the incoming photons that enter and 
are scattered by the local system. 

Each term in ( |24[ ) can be represented diagrammatically. For the case with TV = 5, for example, we plot all different 
classes of diagrams in Fig{^ (We summarize the definitions of all diagrams used in the paper in Figj^. The summation 
in (24) then represents the summation of all such diagrams as shown in Fig{^ each diagram containing a single Green 
function part, and the rest are 5 functions. 


Our main result (24) directly reduces the computation of the iV-photon S matrix to the calculation of the Green 


function of the local system. The result here is related to, but not identical to, the LSZ reduction approach as discussed 
in m- In [12], the iV-photon S matrix is first reduced with the LSZ reduction to the calculation of the Green function 
of the waveguide photons, which is then related to the Green function of the local system after integrating out of 
the waveguide photon fields. In contrast, in our derivation the reduction to the Green function of the local system is 
obtained directly with input-output formalism, which represents a more direct approach. Also, the proof here is very 
general. It uses only the quantum causality relation without any need for knowing the details of the Hamiltonian of 
the local system, which again points to the power of the input-output formalism. 
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FIG. 1: The diagrammatic representation of all six classes of terms that arise in the computation of S matrix for N = 5 
photons. 


V THE CONNECTEDNESS STRUCTURE OF N-PHOTON S MATRIX 


Our main result (21) and (24) reduce the problem of computing fV-photon S matrix, to the calculation of Green 


functions of the local system. Such a calculation, moreover, can be significantly simplified exploiting the connectedness 
structure of the S matrix and the Green function. 

Since the system is time-translation invariant, ,fcjv) must be proportional to 

^ {^^=1 Pi ~ ■ In. general, ,pN',ki,--- ,k]s[) can be expressed as a sum of various terms con¬ 

taining products of several S functions |40j . Among all these terms, we define the term that contains only 
6 (^^^iPi — ^n.d no other 6 functions as the connected Green function, G‘^(pi, • • • ,pjv; ki, - ■ ■ , fc^r). Simi¬ 


larly, as can be seen in (24), the S matrix can also be organized as summing over various terms containing products 


of several 6 functions. Among all these terms, we can again define the term that contains only 6 ~ 

and no other <5 function as the connected part of the S matrix, 


Our main result (24) then immediately implies that for N > \ 




(25) 


and for N = \ 


Spi-ki = 5 {pi - ki) -I- G'^(pi; fci). 


(26) 


Therefore, for A^ > 1, we use the same diagrammatic representation for and G*^ (Figj^. For N = 1, the 
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k 


S(p-k) 



N>1 


K 

^2 

kj^ 


-A 

Pi 

Pn 


P\---Pn'A---K 


■N 


N>1 



PVPN'A---kN 


N>1 


or 

G''iPi---PN,K---k^) 

N>2 


FIG. 2: Various diagrams used in the paper. The corresponding mathematical definition is shown on the right. The legs to 
the left (right) represent input (output) momenta. There is no distinction in the ordering among the input or output momenta. 


diagrammatic representation of is included in Figj^and the result of (26) can then be represented in Figj^ 
Moreover, it is known that a Green function in general can be cluster-decomposed as gOlHI]: 


Mq 

G(pi, • • • ,Pn', ki, - ■■ , Um) = X/X/ iPBi'^ksPi) (27) 

B P i=l 

where ;B is a partition of an ordered list {1,2, • • • , N} into a collection of subsets. (For example, for a list {1, • • • ,5}, 
a partition results in {!}, {2,4}, {3, 5}.) BP is the same partition as B but acts on the permutated ordered list 
{P(l), P(2), • • • ,P{N)}. (For example, the same B in the example above, acting on the permutation {5,3,4, 2,1} 
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FIG. 3: The single-photon S matrix is the sum of a <5 function and a two-point Green function. 


would result in {5}, {32}, {41}.) Mg is the number of subsets and Bi, BPi are the respective i-th subsets. We also 
use the shorthand notations pg. = {pj\j S Bi\ and fcgp. = {kj\j S BPi\. For a given partition, we only sum over 


all the distinct permutations under the exchange symmetries in each sublist. Graphically, (27) states that the Green 


function can be decomposed into a summation of all possible products of the connected Green functions. 


By (25), (26) and (27), we prove in the appendix that our main result (24) implies similar cluster-decomposition 


properties for the S matrix of the system 


Mb 


=EEn ^PBpkBPi ' 

B P i=l 


Below, we give three simplest examples of (28) when N = 1,2,3, respectively: 


(28) 


Sp^-ki — 5{pi — ki) + G{pi \ ki) — . 


(29) 


^piP2',kik2 — ^pi\kiSp2-^k2 4" ^pi',k2^P2',ki + S {pi,p2', fci, fc2) • 


(30) 


5, 


PlP2P3',klk2k3 


'y ^Pi;kp^ Sp2\kp^ ^P3',kp^ 


\_^Pi,P2\kp-^,kp^^P3',kp^ ^P2,P3-,kp2,kp^^Pi-,kp^ Spi,P2\kp-^,kp^^P2\kp^ 


+S^{pi,P2,P3',kl,k2,k3) . 


(31) 


(30) and (31) are represented diagrammatically in Figj^ Using such diagrammatic representation, N photon S matrix 
can be straightforwardly decomposed. The results here thus reduce the computation of the Wphoton S matrix, to 


the evaluation of the connected 2Wpoint Green function of the local system. 
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FIG. 4: The cluster decomposition of two- and three-photon S matrices. 

VI COMPUTE SYSTEM’S GREEN FUNCTION WITH THE EFFECTIVE HAMILTONIAN FOR THE 

SYSTEM 


In this section, we will prove that the Green function of (22) can be computed using the effective Hamiltonian of 
the local system. The main result in this section is 


G{t'i ■ ■■ ■ In) = G{t[ • • • • • • t ^), 


(32) 


where 


G{t\ ■ = (- 7 ) (0|T a{t\) ■ • •a(t(v)a'l'(G) • ■ •at(t7v)|0) 


(33) 


with operators 


d(t) = a , at (t) = at , 


(34) 


where 


• 7 t 

^eff = Hsys ^ 


is the effective Hamiltonian for the system. With the identity (32), the computation of the Green function is simplified 


(35) 


since no operators of waveguide photons are involved. We only need to solve a system which has a finite, and typically 
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small, number of degrees of freedom. 


(32) can be proved in the path integral formulation. The proof is summarized as: 


G{t[ ■ • • • - tjv) 




/ T>[cfc,Cfc,a, 

/ T>[cfe,c^,a, 

J V [a, a*] a{t[) ■ ■ ■ a(t^)a*(ti) • • • a*{tN)e^ 
J T) [a, a*] 


G{t'l ■ ■ ■ t'jy] ti • • • tjv) • 


(36) 


The first step in (36) follows [H] to express Green function (22) by path integral. Here for simplicity, we assume 
the local system is characterized only by a and a*. L is the Lagrangian associated with the full Hamiltonian ([^: 


L = y dk cl{idt - k)ck - ^ J dk {cla + a*Ck) + Lgys, 


(37) 


in which Lgys is the local system’s Lagrangian obtained by Legendre transformation on the system’s Hamiltonian 


77sys ■ 


The second step in (36) is the key of the proof. As was done in the standard approach involving generating 


functional, we introduce the propagator of free waveguide photon 




d^ 


uj -k + i0+ 


(38) 


and integrate out the waveguide degrees of freedom: 


J 'D[ck,cl]e’’-^ y p [cfc, c^] 

= e^fdtL.ys ( 39 ) 


where 


= y dkGf{t-t')=i y y dk 

= TT / — = irSit - t'). 

./ 27r 


— k 


— iTrS{uj — k) 


(40) 


As a result, we obtain the effective Lagrangian 


LefT - Lgys -)- ZTT^ Q, a - Lgyg A ^ 2 ^ ^ * (^^) 

Here, the imaginary part of the effective Lagrangian arises since the waveguide degrees of freedom that we are 
integrating out forms a continuum. 
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In the last step of (36), the path integral to the system’s degrees of freedom with the effective Lagrangian (41) 


corresponds exactly to the alternative Green function (33) in the Heisenberg picture |44] . The effective Hamiltonian 


(35) is obtained from the effective Lagrangian (41) by Legendre transformation. 


Since in practice, all evaluations of the Green function will be carried out using the effective Hamiltonian (35), in 


the following we will no longer make the distinction between the Green function of system defined in the full coupled 
Hamiltonian versus the effective Hamiltonian, i.e. we will no longer make distinction between the left and right hand 


side of (32). In what follows, the time evolution of all system operators is considered as that in (34). 


(32) allows us to compute the full Green function of the local system. To compute the S matrix, we only need the 


connected part of the Green function. In practice, the connected part of the Green function may actually be obtained 
in a simpler fashion without the need to evaluate and perform a cluster decomposition of the full Green function. 
This will be illustrated in the example below. 

VII EXAMPLE: THREE-PHOTON S MATRIX WITH KERR NONLINEAR CAVITY 


As an example of the application of the formalism developed in this paper, we compute the S matrix of three-photon 
transport in a single-mode waveguide side-coupled to a ring resonator incorporating Kerr nonlinear media. The full 
Hamiltonian has the same form as ([^ with the specific form of Hsys'- 

(42) 

where a is the annihilation operator of cavity photon satisfying the standard commutator relation [a, a^] = I. Let 


zj t I A t t 

Hsys = u!c a' a + — a' a'aa 


a = ojc — i'i, the effective Hamiltonian (35) in this case is 


and can be diagonalized as 


Heff = a a^a + ^ a^a^aa , 


Hee\n) = an+^n{n-l) |n) . 


(43) 


(44) 


Use the formalism as described above, we now compute the S matrix for up to three photons. From (29)-(31) we 


only need the connected S matrix for up to three photons, which is the focus of the calculation here. We consider 


single photon S matrix first. Using (32), the two-point Green function is computed as 


G(t';t) = —7^(O|a(t')|n)(n|a^(t)|O)0(t'— t) = —ye — 


(45) 
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where a completeness set of states of the system is inserted. To obtain the final result, we note that only the single 
photon state contributes to the summation. The single photon S matrix is then: 

Sp-k = \ — k) = [1 +Sk]6{p-k), (46) 

fc — Wc + 17/2 

where, for later convenience, we defined 

Sk = -'ijr— ■ ( 47 ) 

k — a 

To compute the connected two-photon S matrix, we first compute the four-point Green function G 2 (t),t 2 ;Gj^ 2 )- 
Depending on the values of the four time labels, the time ordering operation would give rise non-zero terms that can 
be classified into two types: {aa^aa^) and (aaa^a^), i.e. 

2 

(^1 > ^2; > ^2) > (48) 

i=i 

with 

- tpJ0(tp^ -^Q2)^(^Q2 (49) 

P.Q 

G^'^^t[,t2;ti,t2) = (-7)^Xl(®l“(^Qi)“(^Q2)«^(^D)«^(ip2)|0)^(^Qi -^Q2)^(^Q2 -^Pi)^(^Pi -^^ 2 ), (50) 

P,Q 

where both P and Q are permutations over indices {1,2}. We calculate each term by inserting the complete sets of 
eigenstates, which results in: 

(0|a(t})a'l'(ti)a(t2)a'''(t2)|0) = ^ {0\a{t[)\m){m\a\ti)\n){n\a{t2)\l){l\a\t2)\0) 

m,n,l 

= (0|a(t})|l)(l|at(G)|0)(0|a(t')|l)(l|a^(i2)|0) = e—0i-‘i)e-*“02-‘D ^ (51) 

and 


Ht[)a{t'2)aHh)a\t2)\0) = (0|a(t})|l)(l|a(4 )|2)(2|at(G)|l)(l|ani2)|0) 

= 2e“*“0i“‘2)g-h“+x)(*2-*i) 


(52) 


Then by (23), the Fourier transformation of (49) and (50) are 


G^^Hpi,P2]ki,k2) = -^^SpQ^Sk 

^ P.Q 
i 


G^'^\pi.P2\ki,k2) = -^SpQ^Skp^ 


PQ 2 - kp^ - it 

1 


<5(pi +P2 - fcl - ^ 2 ) : 


(53) 


P.Q 


^1 + ^2 - 2a - X 


5 {pi +P 2 -ki- k 2 ), 


(54) 
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where in (531 an infinitesimal imaginary part in the denominator arises due to the Fourier transform of the 9 function. 


Moreover, we note that 


1 


+ iTrS{p — k). 


(55) 


p — k — ie p — k 

On the other hand, since the connected two-photon S matrix contains only a single 5 function, when we apply (55) 


to (531, only the principal part contributes to Therefore, we have 


^?iP2-Mk2 = +P2-ki- fca), 


(56) 


with 


= -—Vs 
P,Q 

<(2) _ * 

1r> • h-. hr, 

TT 


Skf 


. . = - Vsoo Sfe, 

PlP2\kik2 ^ / j PQi 

P,Q 


PQ2 - kP2 
1 


ki + k 2 — 2a — X 

For the connected two-photon S matrix, we can sum over all the permutation terms and obtain a compact form : 


(57) 

(58) 


Sp,p2-Mk2 - (sfci + «fc2) ^ _ 2^ _ ^ 


<5(pi -I-P2 - fci - k2). 


(59) 


The final result (591 indeed has the exact analytical structure constrained by the cluster decomposition principle . 


The only singularities are isolated poles corresponding to one and two-photon excitations in the local system. All 
principle parts cancel. 

Finally, we sketch the computation for the connected three-photon S matrix, we start by computing the six-point 
Green function Similar to the previous calculation on four-point Green function, terms that 

have non-zero contributions to the six-point Green function can be classified into five types: (aa^aa^ao^), (aaa^a^aa^), 
(aa^aaa^a^), {aaa^aa)a))^ and {aaaa^a)a)). Take the last type as an example, we compute its contribution to the 
six-point Green function is 

= (-7)^Xl^®l“(^Qi)“(*Q2)«(^Q3)“^(V)a^(^P2)a^(^P3)|0) 

P,Q 

- ^Q2)^(^Q2 - ^Q3)^(^Q3 “ tPiA^Pi - - V) 

= 6 (— 7 )^ ^ g-*“(‘Qi-*j'3)e“*(“+x)(‘Q2“*^2)e-*(“+2x)(i'Q3-tpi) 

P,Q 

“ ^Q 2 )^(^Q 2 “ ^Q3)^(^Q3 “ ^PAi^Pl - tP2W{iP2 - V) > 


(60) 
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and then its Fourier transformation 


G^^Hpi:P2,P3;ki,k2,k3) = 


^ fci + fc 2 + fca - 3a - 3x Pq^ +PQ^-2a-x kp^ + kp^ - 2a - x 


xS{pi +P2+P3-ki-k2- fca) 

- ^■^‘'pLprMk^kJiPi +P2+P3-ki-k2- ks) 


(61) 


where P, Q are permutations over indices {1,2,3}. Since (611 contains only a single 6 function, all terms in (61) 


contribute to the connected three-photon S matrix. Similarly, we calculate the other four types’ contributions to 


^piP 2 P 3 -kik 2 k 3 by applying (55) and keeping only the principal parts. The final result is summarized as: 


where 


^PlP2P3-Mk2k3 = Hi^plP2P3-Mk2k3^^P^ + P2 + PS ” ^1 - /C2 ” ^ 3 ) , 

i=i 


(62) 


iM 


( 1 ) 


PlP2P3-Mk2k3 4^2 ^kp^+kP3 PQ3 ^>^P3 

P.Q 


(63) 


iM 


( 2 ) 
P1P2P3 


1 ^ VI 

■Mk2k3 - 2^2 _ kp^ p^^ _ 2a- X ’ 


(64) 


. . , ( 3 ) _ 1 

* PlP2P3-Mk2k3 - “^Z^'®P<3i®'=^’2+'=P3-PQ3®'=^’3}};^^7r^^+fcp^ _ 2a - X ’ 
PiQ 


1 


(65) 


•A .<(4) _ ^ ^ ^ 

*-^PiP2P3;fcifc2/=3 - ^2 ^^PQi^k^2+kp3-PQ3^i^^3 p^^ + p^^ _ 2a - X kp^ + kp^ - 2a - X ’ 


( 66 ) 


iM 


(5) Sfo 


P 1 P 2 P 3 


_ 'V' 

•,kik2k3 9 2 Z-^^PQl 


^ ki+k2+k3-3a- 3x PQi -f pq^ - 2a - x kp^ + kp^ -2a- x' 


(67) 


As a check, let x “t oo, only (63) contributes to the connected three-photon S matrix, which agrees with the 


result in the case of single two-level atom m, as expected. Also, according to the cluster decomposition principle, 
the connected three-photon S matrix should only contain poles corresponding to the single, two- and three-photon 


excitations [45]. In (63l-(67), in addition to various poles corresponds to single, two- and three-photon excitations, 


there are also various singularities associated with the principal parts. One can actually prove that these principal parts 
cancel each other when all permutations are summed together. A systematic treatment of such analytic properties of 
the connected A^-photon S matrix is beyond the scope of this paper and will be carried out in future works. 
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SUMMARY AND FINAL REMARKS 

To summarize, in this paper, using the input-output formalism, we provide a computation of 7V-photon S matrix 
in waveguide QED systems. The main result here is the connection between the A-photon S matrix and the Green 
function of the local system. We also discuss the connectedness structure of the S matrix and the Green function, and 
how such structure can be used to simplify the computation. Our results are applicable independent of the details 
of the local system’s Hamiltonian, and therefore point to some universal aspects of the properties of waveguide QED 
systems. As a computational tool, the results here lead to a powerful scheme for A-photon S matrix calculation. 

Aside from waveguide QED systems, understanding the scattering property of a local quantum system coupled to 
a continuum has been a problem of fundamental importance in many other branches of physics. For example, in 
condensed matter physics, the transport properties of a quantum dot can also be formulated in a similar fashion [46) . 
We therefore expect our development here to be useful beyond waveguide QED systems. 


ACKNOWLEDGEMENT 


This research is supported by an AFOSR-MURI program. Grant No. FA9550-12-1-0488. 


APPENDIX 


Our aim here is to prove the cluster decomposition property of the A-photon S matrix (28). 


We start from (24), and use the cluster decomposition property of the Green function (27) to expand every term 


in the summation in terms of sum over products of connected Green functions and ’’bare” 6 functions (i.e. the S 
functions that arise explicitly in ( |24[ )). Among all these terms resulting from the expansion, we consider the term as 
represented by the diagram in the Fig[^ (a). This diagram contains a sub-piece with M legs. Within the sub-piece 


all connected parts are 2n-point Green functions with n > 2. By (25) this sub-piece is already in the form of the 


products of connected S matrices. The rest of this term contains only bare 5 functions. This term obviously arises 
from the expansion of a term in ( |24[ ) containing a 2M-point Green function. 

We combine the term as shown in Fig[^ (a), with the terms shown in Figj^ These terms were chosen from the 


expansion of the terms in (24) containing 2M'-point Green function with M' > M. In selecting these terms, we keep 


the M-leg sub-piece identical in all the terms, and choose only terms where the remaining sub-pieces are either bare 
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FIG. 5: (a) A term in (241 with a M-leg sub-piece that contains only the product of the connected 2n-point Green function 

with n > 2. (b) Same as in (a), but with all the 5 functions replaced by single photon S matrices. This diagram now contains 
the product of the connected S matrix. 


6 functions or two-point Green functions. Repeatedly using (26), the summation of all these terms then results in a 


diagram that has the same structure as Figj^ (a), but with the bare 6 functions all replaced by the single-photon S 
matrices, as shown in Figj^b). The summation therefore results in a product of the connected S matrices. 

Since after the expansion of ( [24| in terms of the product of the connected Green’s function, every term shows up 
once and only once in a summation of the form shown in Figj^ (a) and Figj^ We have therefore proved that (24) can 
be summed to give the cluster decomposition property of 
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